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ABSTRACT 

In this work, we investigate the nature of the host galaxies of long Gamma- Ray bursts 
(LGRBs) using a galaxy catalogue constructed from the Millennium Simulation. We 
developed a LGRB synthetic model based on the hypothesis that these events originate 
at the end of the life of massive stars following the collapsar model, with the possi- 
bility of including a constraint on the metallicity of the progenitor star. A complete 
observability pipeline was designed to calculate a probability estimation for a galaxy 
to be observationally identified as a host for LGRBs detected by present observational 
facililties. This new tool allows us to build an observable host galaxy catalogue which 
is required to reproduce the current stellar mass distribution of observed hosts. This 
observability pipeline predicts that the minimum mass for the progenitor stars should 
be ^ 75 M0 in order to be able to reproduce BATSE observations. Systems in our 
observable catalogue are able to reproduce the observed properties of host galaxies, 
namely stellar masses, colours, luminosity, star formation activity and metallicities as 
a function of redshift. At 2; > 2, our model predicts that the observable host galaxies 
would be very similar to the global galaxy population. We found that ~ 88 per cent of 
the observable host galaxies with mean gas metallicity lower than 0.6 Z0 have stellar 
masses in the range lO^'^-lO^'^'^M© in excellent agreement with observations. Inter- 
estingly, in our model observable host galaxies remain mainly within this mass range 
regardless of redshift, since lower stellar mass systems would have a low probability 
of being observed while more massive ones would be too metal-rich. Observable host 
galaxies are predicted to preferentially inhabit dark matter haloes in the range 10^^- 
IO^^-^Mq, with a weak dependence on redshift. They are also found to preferentially 
map different density environments at different stages of evolution of the Universe. 
At high redshifts, the observable host galaxies are predicted to be located in similar 
environments as the global galaxy population but to have a slightly higher probability 
to have a close companion. 

Key words: gamma-rays: bursts - methods: numerical - stars: formation - galaxies: 
evolution, interactions. 



1 INTRODUCTION 

Gamma-ray bursts (GRBs) are brief pulses of 7-ray radia- 
tion observed on average once a day at random directions in 
the sky. They are the brightest sources in this region of the 
electromagnetic spectrum, and they have been systemati- 
cally studied in the past two decades (e.g. iMeszaroa |2006| . 
and references therein). The origin of GRBs is cosmologi- 
cal, as determined from the measurement of their redshifts 
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(e.g. iMetzger et alj 



iBloom. Diorgovski &: Kulkarnil 

I2OOII : iBloom et al.l [2003J), and in ma ny cases, their 
host galaxies have been identified (se e Le Floc'h et al.l 
I2OO3I : ISavaglio. Glazebrook fc Le Borgnj 12003 . and refer- 
ences therein). Their distances imply the release of large 
amounts of energy (~ 10^^ ergs) in a short time-scale, which 
suggests that the origin of these phenomena could be as- 
sociated to_the accretion of matte r onto a compact object 
llWooslevI I1993I: iFrver et all Il999l: iMacFadyen fc WooslevI 



ll999l : |Panaitescu fc Kumaill200ll : lFrail et al.lboOll ') 



Two populations of G RBs are apparent from the dis- 
tribution of their duration (|Kouveliotou et al.lll993l ). Those 
lasting less than 2 s are known as short GRBs, while 
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longer events are called long GRBs (LGRBs). These are 
more frequently observed and more precisely located, which 
makes their properties better known. LGRBs are always 
fou nd in host galaxies with ongoing star formation activ- 
itv llLe Floc'h et al.ll2003l:IChristensen. Hiorth fc Gorosabell 



I2OO4I : IProchaska et al.M2004h . and some of them were ob- 
served to be associated with Type Ib/c core-coll apse su- 
pernovae taala ma etall Il998l : iHiorth et al.1 1200^1 . These 
observations suggest that the progenitors of LGRBs are 
massive stars. Stellar evolution models have tried to pro- 
vide a consistent scenario where a LGRB can develop. 
The result is the so called collapsa r model (|Wooslevl 
ll993l : lFrver. Wooslev fc Hartmannll 19991 ') . in which this phe- 
nomenon is produced during the collapse of a massive star, 
due to the accretion of part of the envelope onto the recently 
formed black hole. Nevertheless, the properties of LGRB 
progenitors (mass, metallicity, rotation velocity, binarity, 
etc.) are still a matter of discus sion. LGRBs h ave been as- 
sociated to single massive stars JWooslev 19931) . perhaps in 



meta llicity biased environments IjHirschi. Mevnet fc Maedei 



20071 ). 



2005 ). and also to stars in binary systems (jBissaldi et al 



The importance of understanding the nature of the 
LGRB progenitors is beyond the interest of only stellar 
evolution and black hole formation. Given their connec- 
tion to massive stars and their large luminosities, LGRBs 
might be a powerful tool to investigate the star forma- 
tion in the early Universe, at r edshifts for which stan- 
dard tools become ineffective (e.g. IWiiers et al.lll998l ). The 
knowledge of the properties of the stellar progenitors would 
allow the assessment of possible biases originated when 
LGRBs are used as tracers of star formation. Several au- 
thors suggest that the cosmic LGRB rate does not follow 
the s tar formation rate rPai me. Rossi fc Mochkovitchll200q : 
ISalva terra fc Chincarini 20071) ■ The interpretation of their 
data requires the assumption of a differential evolution of 
the comoving LGRB rate density with respect to the co- 
moving star formation rate ( SFR) density. The same conclu- 
sion arises from the work of IWolf fc Podsiadlowskil (|2007|), 
who show that the luminosity function of LGRB hosts dif- 
fers from that of core-collapse supernovae hosts, which are 
considered unbiased tracers of star formation. According to 
these authors, the LGRB host luminosity function can be 
reproduced by requiring LGRB progenitors to have metal 
abundances lower than that of the Sun. 

An indirect procedure to investigate LGRB progeni- 
tors is to characterise the stellar populations of their host 
galaxies. Given the connection between LGRBs and star 
formation, galaxy formation models can be used to assess 
the validity of LGRB progenitor models, comparing their 
predictions about the properties of these stellar popula- 
tions to host gal axies ob s ervatio ns. Recent works such as 
that of iKauffma nn et al.l l|2004) . suggest that there is a 
strong relation between star formation and nuclear activity 
in a galaxy with its environment, at constant stellar mass. 
Hence, the environment of host galaxies might constitute 
an independent probe for progenitor models. Observational 
works on host galaxies environment are still inconclusive on 
wheth er host galaxies inhabit regions with certain character- 
istics (|Fvnbo et al.ll2002 ':ljakobs son et aLll2005l : lThone et all 
120081 ). iBornancini et al. (2004) s uggest that host galaxies are 
field galaxies, while lWainwright. Berger fc Penprasd (|2007l ) 



find that a considerable fraction of their sample of 42 host 
galaxies shows evidence of interaction with other galaxies. 
The question of whether there is a connection between the 
occurrence of a LGRB and the local density of galaxies re- 
mains unanswered. 

Since galaxy formation is a highly non-linear process, 
the properties of host galaxies predicted by differ- 
ent progenitor scenarios are better studi ed by means 
of c os mological numeri ca l sim ul ations (jKatz fc GunnI 



I991I: 



iNavarro fc White! 
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19931: iMosconi et a" 
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Springel fc Hernquist 2003; Scannapieco et al 
20061 ). ICourty, Biornsson fc Gudmundsson (2004) 



2001 



20051 . 



hydrodynamical simulations of structure formation to iden- 
tify a galaxy population whose properties reproduce the 
observed ones. According to these authors, the luminosity 
distribution of host galaxies is reproduced if galaxies are re- 
quired to have high star formation efficiency. However, their 
simulations do not make a ny prediction on t he properties 
of the LGRB progenitors. iNuza et al.l (J2007|) developed a 
Monte Carlo code to simulate the production of LGRBs in 
hydrodynamical simulations of galaxy formation, assuming 
that their progenitors are massive stars as proposed by the 
collapsar model. These authors follow the LGRB produc- 
tion as the structure forms and evolves in the Universe. 
Their results suggest that LGRB progenitors would be 
low metallicity stars {Z < 0.3Zq), and hence LGRBs 
would be biased tracers of star formation principally at low 
redshift. However, their simulations explore a small volume 
of (10h~^ Mp c )*^ do minated by field galaxies. Recently, 
ICampisi et al.l (|2009l ) have constructed a simulated popula- 
tion of host galaxies based on the semi - analyt ic model of 
galaxy formation of iDe Lucia fc BlaizotI (120071) applied t o 
cosmological simulations developed by IWang et al.l l|200a ) . 
These simulations describe a larger volume of the Universe 
a square box of 125h~^ Mpc in a side). ICampisi et al.l 



2009D explore three LGRB progenitor models based on the 
collapsar model, one of them assuming a given stellar mass 
value and an age cut-off and the other two also including 
metallicity cut-offs. These authors find that models with 
very low metallicity progenitors {Z < O.IZq) could explain 
the luminosities, colours and metallicities of the observed 
host galaxies. Their results support previous claims of 
LGRBs being biased tracers of the star formation. 

Clearly, the use of larger simulations with better resolu- 
tion and more detailed descriptions of the physical processes 
driving the dynamical, chemical and star formation histories 
of galaxies allows a better modelling of the host galaxies of 
LGRBs. However, an important aspect that should not be 
disregarded is that, for a proper comparison with observed 
host galaxies, the effects of the detectability of these galax- 
ies must be taken into account. And this is the main con- 
tribution of our work where we developed an observational 
pipeline which allows us to mimic, at least in part, biases 
affecting the observations. With this new tool we build an 
observable host galaxy catalogue which can be more fairly 
compared to current available observations. 

In this work, we develop a semi-analytical model for the 
host galaxies of LGRBs adopting the collapsar model for 
the progenitor stars. We apply it to one of the largest cos- 
mological simu l ations available, the Millennium Simulation 
ISpringel et al.l ((20051). The Millennium Simulation follows 
the evolution of dark matter in a box of 500h~^ Mpc in a 
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side, hence providing a cosmolog ically representativ e volume 
a factor of 64 larger than that of Wang et alj (120081 ) and fac- 
tor of ~ 10^ larger than that of^Nuza et al.' ('2007'). Wc take 
the galaxy catalogue built up by De Lucia fc Blaizot ( 2007), 
who describe the evolution of baryonic matter including star 
formation, active galactic nuclei and supernova feedback. 
Using the star formation and chemical properties of galax- 
ies in this catalogue to implement different LGRB progeni- 
tor scenarios, we determine the properties of the overall host 
galaxy population and its environment. We compute the de- 
tectabiUty of the host galaxies in order to compare the pre- 
dictions of these sce narios to the observations compiled by 
ISavaglio et~all (|2009l ). 

This work is organized as following. In Section [2l 
we briefly describe the cosmological simulation and semi- 
analytic galaxy catalogue used in our work. In Section [31 we 
develop our implementation of LGRB progenitor scenarios. 
Sections [4] and [5] show our results on host galaxy properties 
and environment, respectively. Finally, Section [S] presents 
our conclusions. 



2 MILLENNIUM SIMULATION 



The Millennium Simulation (jSpringel et al.ll2005h is one of 
the largest cosmological simulations of structure formation 
publicly available. It describes the evolution of the cold 
dark matter while the baryonic matter can be included 
via semi-analytic models on top of the numerical simula- 
tion. The Millennium Simulation describes the formation 
of structure tracking ~ 10^° dark matter particles of mass 
8.6x10^ h~^ M0 distributed in a cubic region of 500 h~^ Mpc 
side, using the ACDM cosmogony. The adopted cosmolog- 
ical parameters are Qm ~ 0.25, f2b ~ 0.045, f^A = 0.75, 
as = 0.9 and n — 1, where the Hubble constant is Hq = 
100/ikms~^ Mpc~^ with h — 0.73. These parameters are 
chosen in consistency with results obtained by the joint 
analysis of the Two-degree Field Galaxy Redshift Survey 
(2dFGRqj) and the Wilkinson Micro wave Anisotropy Probe 
(WMAP) data (jSanchez et al.ll2006l ). 

The simulation was performed using a modified ver - 
sion of GADGET-2 hydrodynamical code (|Springell [20051 ). 
A friends-of-friends (FOF) algorithm identified non-linear 
dark matter haloes within the simulation in an automatic 
manner. Two particles belong to the same halo if their sep- 
aration is less than 0.2 times the mean separation between 
particles. Only groups of more than 20 particles are identi- 
fied as haloes. This restriction sets a lower limit on the mass 
of a halo of 1.72 x 10^" /i~^ M©. Substructures orbiting within 
each virializcd ha lo are ident ified applying a SUBFIND al- 
gorithm (Springcl et al.l 120011 ). Each dark matter halo has 
a central (Type 0) gala:xy and one or more satellite galax- 
ies. Satellites were at some point central galaxies of smaller 
haloes which suffered a merger with the halo their currently 
inhabit. There are two types of satellite galaxies, Type 1 
galaxies are located at the centre of a subhalo associated 
with a FOF group, while Type 2 galaxies have lost their 
dark matter subhalo after falling onto a more massive halo. 



http://www.inso.anu.edu.au/2dFGRS/ 



Mean properties of synthetic galaxies (e.g. SFR, stel- 
lar mass, gas mass, metallicity) are obtained by applying 
semi-analytic models t o the struc t ure fo rmation simulation 
(White fc Frcnk lOOll: ICold [l99ll : iLacev fc Silk 1991 ). Es- 
sentially, semi-analytical models describe the collapse of 
baryonic matter following dark matter haloes. After col- 
lapsing, the gas cools and gravitational instabilities produce 
episodes of star formation. To reproduce observations, mod- 
els incorporate photoionization processes in the intergalac- 
tic medium, the growth of supermassive black holes during 
galactic mergers, supernova and AGN feedback, star forma- 
tion rate enhancement in mergers, and formation of heavy 
elements for each formed stellar populat ion. The galaxy cat- 
alogue we use in this work is that of lOe Lucia &: Blaizotl 
(|2007h . 

There is recent evidence for an excessiv e red- 
dening of Type 2 galaxi es (jWeinmann et al.l 120061 : 
IPerez. Tissera fc Blaizoj |2009| ) which might be due to the 



poor physical treatment of this type of galaxies. In order 
to avoid spurious trends, we exclude Type 2 galaxies from 
our analysis in Section and Section [J] However, Type 2 
galaxies ought to be included in the environmental analy- 
sis performed in Section [5] in order to correctly trace the 

underlying mass distribution^ 

Hence, the catalogue of lDe Lucia fc Blaizoj (|2007l ) pro- 
vides us with the spatial galaxy distribution, their stellar 
masses, dark matter haloes, star formation activity, colours, 
luminosities and mean metallicities. Regarding the latter, 
the public catalogue makes available the mean metallic- 
ity of the cold gas component and of the stellar popula- 
tion as a whole (i.e. averaged over new and old stars in 
a galaxy). Therefore, as a proxy for the metallicity of the 
LGRB progenitors, we take the mean metallicity of the cold 
gas component at the time th ese stars were born. Because 
the lDe Lucia fc Blaizod (|2007|) model assumes instantaneous 
recycling, this metallicity is slightly higher than that of 
the newly born stars, but this correction is negligible com- 
pared to the uncerta inties in the measured metallicities (e.g. 
ICampisi et~al]|2009l ). 



3 SCENARIOS FOR LGRB PROGENITORS 

3.1 Intrinsic LGRB rate 

We consider two scenarios for LGRB pr oduction , both 
of them based on the col l apsar model (IWooslevI Il993l : 
iMacFadven fc Wooslevlll999l : iFrver et al.lll999l ). In our sce- 
nario I we take as LGRB progenitors all stars above a certain 
minimum mass nimin, with no other restriction whatsoever. 
In this scenario LGRBs are unbiased tracers of star forma- 
tion. Accordingly, we obtain the intrinsic LGRB rate in a 
given galaxy g at a particular redshift z as 



rGRB(ff,z) = SFR(5,2) 



/, 



100 M0 



£,{m) dm 



L .. mt(ra) dm 



(1) 



where 0.1 M© and 100 Mq are the lower and upper mass cut- 
oflts of the IMF ^(m) given bvlchabricr (2003). According to 
Eqn.[T] the production of LGRBs is not delayed with respect 
to the starburst that created the progenitor stars, which is 
justified because in the collapsar model only massive stars 
are LGRB progenitors. 
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In our scenario II, we assume that only massive stars 
(m > rrimin) below some metallicity threshold (Zc) are 
able to prod uce LGR Bs, as in the progenitor models of 
iHirschi et all (20g5_) or JYoon. Langer fc NormanI (|2006). As 
explained before, we take the mean metallicity of cold gas 
in each galaxy as representative of the metallicity of the 
progenitor stars. The consequences of this hypothesis will 
be discussed in later sections. In this scenario, the intrinsic 
LGRB rate is given by Eqn. [T] for galaxies g with metal- 
licity Zg < Zc, and is null for those with Zg > Zq. 
Three realizations of this scenario were computed, adopting 
Zc — 0.1, 0.3, 0.6 Z0 (scenarios II. 1, II. 2 and II. 3, respec- 
tively). 

3.2 Observed LGRB rate 

Given a metallicity threshold, we adjusted the only free pa- 
rameter of our scenarios (mmin) to reproduce the LGRB rate 
measured by the BATSE exp eriment. For this p urpose, we 
took the off-hne GRB search of lStern et all l|20011 ). which de- 
tected 3475 GRBs with Tgo > 2 s during a live-time of 6.37 yr 
(70 per cent of the 9.1 yr that BATSE was active), scanning 
67 per cent of the sky. This translates into a full-sky LGRB 
rate of 814 yr~^ above the threshold of the off-line search. 
The choice of this particular experiment was motivated by 
its good statistics and the availability of an accurate model 
for its detection efficiency. 

To obtain the observable rate predicted by a given sce- 
nario, we first compute the comoving LGRB rate density for 
each redshift as 



nGRB(2:) = V ^ y^rGRB(ff,^), 



(2) 



where V = (500 hT^ Mpc)^ is the comoving volume of the 
Millennium Simulation. The full-sky observed LGRB rate is 
then 



^GRB = 







nG-RS.{z) dV 

—77—. ^Pdct.BATSEi^J -7- , 

(1 -I- 2) dz 



(3) 



4-Kcdf{z) 



(4) 



where 

H ^ H{z){i + zy 

is the derivative of the comoving volume with respect to z at 
fixed solid angle, Zmax is the maximum redshift of the simula- 
tion, dL the luminosity distance and H the Hubble constant 
at redshift z, and pdet.BAxsE the probability of detecting a 
LGRB at redshift z with BATSE (i.e. the probability that 
a particular LGRB has a peak flux above the experiment 
threshold). The integration was performed numerically, and 
the described procedure iterated over mmin, until agreement 
with the observed rate was attained. 

The value of Pdct,BATSE(z) depends on the LGRB lu- 
minosity function and spectrum. It was computed using a 
Monte-Carlo scheme to simulate, at a given redshift, a large 
population of LGRBs with different luminosities and spec- 
tra. Then, their photon fluxes in the BATSE energy band 
were estimated. For each LGRB in this population, a second 
Monte-Carlo procedure rejected those events which would 
be un detectable, taking into account the trigger efficiency 
of the lStern et all (|200ll ') off-line search. If the LGRB emis- 
sion were isotropic, the fraction of retained LGRBs would di- 



Table 1. Main characteristics of our scenarios for the simulated 
host galaxies of LGRBs. Column (1) gives the name of the sce- 
nario. Column (2) lists the minimum stellar mass obtained for 
progenito r stars of LGR Bs to reproduce BATSE observations, 
adopting IChabrieJ 1I2OO3I ') IMF. Column (3) gives the cold gas 
maximum metallicity cut-off. Column (4) shows ratio between 
the percentage of simulated host galaxies over the mass range 
fa io8-5-l 0.3jy^^ ^j^j ^^g corre sponding value obtained from the 
sample of lSavaglio et al.l 1120091 ). 



Scenarios mmin (M0) Zq Rm 



I 


91.4 ±0.1 


- 


0.82 


1.1 


13.9 ±0.2 


0.1 


0.02 


11.2 


44.3 ±0.4 


0.3 


0.51 


11.3 


76.0 ±0.3 


0.6 


1.03 



rectly give the detection probability due to the off-line search 
thre shold. To account for the beamed emission of LGRBs 
fe.g. lYonetoku et al.ll2005l ). and assuming that the distribu- 
tion of jet opening angles is independent of the LGRB lumi- 
nosity, spectrum and redshift, we obtain Pdct.BATSE by mul- 
tiplying this fraction by the mean beaming fraction of the 
jets pjet ~ 10"'^. The luminosity func tion and spectral pa - 
rameters distributions were taken from lDaigne et al.l (|2006r ). 
This procedure not only allows us to mimic the obser- 
vational process but also has provided us with estimates of 
the minimum mass for the progenitor stars. The values of 
rrimin obtained, in solar masses, are 91.4 ± 0.1, 13.9 ± 0.2, 
44.3±0.4, and 76.0±0.3 for scenarios I, II.l, II.2 and II.3, re- 
spectively; the uncertainties reflect Poissonian errors in the 
number of observed LGRBs (see also Table 1). 



3.3 Host galaxies 

A meaningful definition of a LGRB host galaxy in our sce- 
narios is not as straightforward as it might seem at first 
sight. The naive definition of a host galaxy, at a given 
redshift z, as being any galaxy g with a LGBR rate of 
''grb(<7,z) > (or identically, SFR((;, z) > 0) is not useful 
because the resulting host population would include galax- 
ies with arbitrarily low LGRB rates. Low rates are better 
understood in statistical terms, as very low probabilities of 
producing a LGRB per unit time. This implies that the cor- 
responding galaxies have a low probability of being observed 
as host galaxies. Hence, this naive definition would generate 
a host galaxies sample biased to low SFR galaxies. To fix 
this problem, the definition of a host galaxy could be based 
on the number of LGRBs N{g, z) = raKsig, z)At produced 
in each galaxy during a time interval At in its rest frame, 
so that host galaxies are only those with N{g,z) > 1, as 
in Campisi et al. (2009). However this cut-off, and the re- 
sulting population, would be dependent on the more or less 
arbitrary choice of At. 

Given the above arguments, we preferred instead a 
probabilistic approach, defining the likelihood of a galaxy 
being observed as a LGRB host. This likelihood is then used 
as a weight to compute the properties of the observable host 
galaxies sample, which would be in this way comparable to 
the observed sample. For a galaxy to be detected as a host at 
least one LGRB must be detected within it by a high-energy 
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observatory. Then, the galaxy itself must be detected, usu- 
ally by optical/NIR telescopes. We would treat the biases 
introduced by each of these observations separately. 

To model the first bias, we compute the probability of 
detecting at least one LGRB within it by any of the high- 
energy observatories monitoring GRBs. Given that event de- 
tection processes follow Poissonian statistics, if RGRB,i{g, z) 
is the contribution of galaxy g at redshift z to the LGRB 
rate observed by experiment i, the probability of observing 
at least one LGRB in a galaxy is 



PHG(g,z) = 1 - exp 



-^-Rgf 



^AT^ 



(5) 



where ATi is the time interval of the LGRB search con- 
ducted by experiment i. Eqn. [5] assumes that the searches 
by different observatories are independent. Previous works 
estimate the intrinsic LGRB rate in a t ypical galaxy to 
be of the order of 10~^yr~^ JFrver et al.l il999'). hence the 
LGRB rate observed in the galaxy by any mission must 
be lower than this value. Given that current searches for 
LGRBs have spanned a few years, ^^ 7?GRB,iAT'i <C 1, and 
Phg(<7, -z) — X/ ^GRB.i ATi, which means that the likelihood 
of a galaxy being observed as a host is proportional to its 
contribution to the mean observed number of LGRBs, 



, s rGRB{g,z) dV ir-^ AQi 



(6) 



where An^ is the sky coverage of the experiment i. Note 
that, for a fixed z, cosmological and observatory dependent 
factors in Eqn. [5] become constant, independently of the 
number of observatories considered, and cancel out in the 
normalization of the weights. This means that when com- 
puting mean HG properties such as mass or SFR as a func- 
tion of z, the first bias can be modeled simply by weighting 
the corresponding properties of each galaxy by its intrinsic 
LGRB rate. 

In the computation of the integrated properties of the 
whole population of host galaxies (like the integrated mass 
distribution, for example), also the effects of volume vari- 
ation, time dilation and detectability of the different ob- 
servatories must be taken into account (Pellizza et al. in 
preparation). In this case it is better to compute first the 
mass distribution observed by each experiment at each red- 
shift z, weighting the galaxies in the corresponding snap- 
shot by pHGig,z). Second, the integration in z can be per- 
formed for each experiment, and the resulting distribution 
can be normalized. This has the advantage of avoiding the 
use of the values of ATi and A^li , which are poorly known 
for some of t he observatorie s that detected the LGRBs in 
the sample of ISavaelio et al.l l|2009r i. Only pdct,i{z) for each 
experiment is needed, which is computed as described in 
Sect. \3l2\ for Pdet,BATSE(z)- Finally, these distributions can 
be combined into a single one by adding them, previously 
scaled to the number of LGRBs detected by each observa- 
tory. The r elevant data f or com puting the d etectability were 
taken from Istern e t al.' ('20 0l|) for BATS E, iGuetta fc PiranI 
ll2007fl for Swift Fro nter a et all (120091') for Beppo-SAX 
iPelangeon et alj (|2008l ) for HETE-2 a.nd \Um\ev et all (|l992l ') 
for Ulysses. Thes e experiments d e tected 36 of the 38 LGRBs 
in the sample of ISavaglio et al.l ([20091). For statistical rea- 
sons, only experiments that detected at least five LGRBs in 



the sample were considered. Konus- Wind and NEAR-XGKS 
were discarded because all their LGRBs were observed also 
by other experiments above. To be consistent, the observed 
mass distribution to which the model predictions were com - 
pared was constructed from the data of lSavaglio et al.l (|2009l ) 
by adding the individual distributions observed by each ex- 
periment. 

The second bias is more difficult to model. Given that 
the search for host galaxies is guided by the discovery of the 
LGRBs themselves and usually done with a variety of differ- 
ent telescopes and detectors, in different bands of the elec- 
tromagnetic spectrum and with different sensitivities, the bi- 
ases introduced are unclear. Galaxies with surface brightness 
below the detectability threshold of available instruments 
could introduce a bias towards low metallicity galaxies. An 
extra problem could be caused by dust obscuration which 
could affect the detection of the afterg lows and produce a 
bias towards high metallicity galaxies (|Fvnbo et al.ll2009l ). 
Hence, at least these two effects might combine themselves 
to determine the detectability of a host galaxy. Considering 
also the fact that we have only access to the public galaxy 
catalogue of the Millennium Simulation which provides only 
mean global properties and magnitudes, we adopt the ob- 
served integrated stellar mass distribution as a tool to apply 
the combined effect of observational biases to the simulated 
galaxy sample. We use the fact that the observed integrated 
stellar mass distribution has been affected by observational 
biases although we cannot disentangle their individual ef- 
fects. Hence, we require the simulated LGRB hosts to repro- 
duce the observed stellar mass distribution in order to build 
up the observable simulated LGRB hosts. The procedure is 
explained in detail in next section. Note that, hereafterin, 
we will discuss the trends of the observable simulated hosts 
and the general galaxy population. The former can be com- 
pared to the current observed hosts, but, if this later sample 
changes due to better or different observational techniques, 
our simulated sample should be also consistently modified 
to match the new observed stellar mass distribution. 



4 LGRB HOST PROPERTIES 

The largest and most comprehensively uniform sample of 
host gal ajcy properties availa ble at present is that com- 
piled bv lSavaglio et al.l (|2009l ln. Stellar masses, star forma- 
tion rates, metallicities, absolute magnitudes and colours of 
46 observed host galaxies up to z ~ 3 were obtained by 
these authors comparing their spectral energy distributions 
to those of synthetic stellar populat ions, using the method 
described in iGlazebrook et all (|2004l '). 

A key point of our models is that we require the stellar 
mass distribution of the predicted host galaxies to match 
that of the observed host galaxies. The stellar mass is 
adopted as the property to be reproduced by the mod- 
els since it is now widely accepted that stellar mass is 
a m ore fundamental quant ity for galaxies than luminosity 
(e.g. JKauffmann et al.ll2004l ). It could be possible that those 
LGRB events with no detected hosts occured in very low sur- 
face brightness galaxies and hence, with low stellar masses. 



Data available at http://www.grbhosts.lnfo| 
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Fi gure 1. Stellar mass d istribution of host galaxies in the sample 
of ISavaglio et alj 1 12009) (shaded in gray) , together with those 
predicted by scenarios I (dashed-triple dotted line), II. 1 (dotted 
line), II. 2 (dashed-dotted line) and II. 3 (solid line). 



As mentioned before the effects of dust could also prevent 
the detection of events in high metallicity (and dusty) galax- 
ies. In our models these host galaxies exist (and are part of 
the global galaxy population) but the observability cut-off 
has been determined by using the current observed stellar 
mass distribution. Fig. [T] shows the distribution of stellar 
m asses of observed host galaxies constructed from the data 
of ISavaglio et al.l ([20091), together with those predicted by 
our scenarios. 

As it can be seen from Fig. [T] the scenario which best 
reproduces the observed stellar mass distribution is that 
with Zc = 0.6 Z© (scenario II. 3, see also Table 1). Lower 
metallicity thresholds predict lower observed stellar masses 
for the host galaxies, while no metallicity threshold (sce- 
nario I) predicts larger ones. In the sample of lSavaglio et al.l 
l|200g ). we find that 85 per cent of the studied galaxies 



have stellar masses over t he range ~ 10 



Mq (see also 



ICastro Ceron et al.ll2008l ). The probability of getting a host 
galaxy within this stellar mass range in scenario II. 3 is 88 
per cent, while for scenarios I, II. 1 and II. 2 it is 70, 2, and 43 
per cent, respectively. Then we conclude that scenario II. 3 
predictions agree fairly well with observations while others 
fail, and therefore in the following sections, we will focus 
only on this scenario. 

We stress the fact that its predictions include the eflects 
of host galaxies observability, hence, a proper comparison 
with observations can be made. In order to contribute to 
the understanding of the nature of LGRB host galaxies, we 
will also compare these predictions with the properties (not 
weighted by host galaxy observability) of both the sample 
of all galaxies with mean cold gas metallicities below 0.6 Z© 
(hereafter low metallicity sam ple) and the complete gala xy 
population of the catalogue of iDe Lucia fc Blaizoti l|2007h . 



4.1 Stellar mass 

As a first step towards understanding the nature of host 
galaxies, we analyse their stellar masses as a function of red- 
shift. As shown in Fig. [5J the mean stellar mass of the host 
galaxies as a function of redshift predicted by scenario II. 3 
reproduces the observed mean trend quite well. From this 



Figure 2. Mean stellar mass of the host galaxies as a function of 
redshift predicted by scenario II. 3 (dotted line) with its standard 
deviation (shaded grey band). Open d i amon ds represent the ob- 
served host galaxies of ISavaelio et al.l 112009) . while filled circles 
correspond to their mean values in redshift intervals of 0.5. For 
comparison, we also include the corresponding mean trends for 
the low metallicty sample (dashed line) and the complete galaxy 
population (solid line). 



figure, we can also see that host galaxies are, on average, 
more massive than galaxies in the complete galaxy popula- 
tion while the latter are more massive than those in the low 
metallicity sample. This can be understood taking into ac- 
count that the host galaxy observability is a strong function 
of the star formation rate, an d that the complete galaxy 
population in the catalogue of iDe Lucia fc Blaizot (12007 ) 
follow s a mass-metallicity relationship (e.g. IPe Rossi et al.l 
120091 ). Then, the cut-off adopted for the mean cold gas 
metallicity to reproduce the observed stellar mass distribu- 
tion implies a cut-off in stellar mass since low metallicity 
galaxies are, on average, less massive than the general galaxy 
population (see also Fig.[TJ. However, as the observability of 
a host galaxy depends strongly on its star formation activ- 
ity and most of the small galaxies have low star formation 
rates, the observable host galaxies tend to be, on average, 
the more massive ones among them. As a result, our observ- 
able sample tends to be populated by systems more massive 
than those in the low metallicity sample or in the complete 
galaxy population. 



4.2 Star formation rate 

In Fig. [3] we display the mean SFR of host galaxies as a 
function of redshift predicted by scenario II. 3, together with 
the corresponding mean values for the low metallicity sam- 
ple and the complete galaxy population. As it can be seen, 
the prediction of scenario II. 3 reproduces very well the be- 
haviour of the observed host galaxies. These have higher 
SFR than the mean of the complete galaxy population, and 
much higher than that of the low metallicity sample. The 
good agreement between our scenario II. 3 and observations 
suggests that the observed host galaxies are biased towards 
galaxies with stellar masses in the range 1O^~^°M0, high 
star formation activity and relatively low gas metal content, 
compared to the mean properties of the complete galaxy 
population at a given z. 
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Figure 3. Mean SFR of the host galaxies as a function of redshift 
predicted by scenario II. 3 (dotted hne) with its standard deviation 
(shaded gr ey band). Open diam onds represent the observed host 
galaxies of lSavaelio et al.l 1120091 ). while filled circles correspond to 
their mean values in redshift intervals of 0.5. For comparison, we 
also include the corresponding mean trends for the low metallicty 
sample (dashed line) and the complete galaxy population (solid 
line) . 



Figure 4. Mean SSFR of the host galaxies as a function of red- 
shift predicted by scenario II. 3 (dotted line) with its standard 
deviation (shaded grey band). Open d i amon ds represent the ob- 
served host galaxies of ISavaelio et al.l 112003) . while filled circles 
correspond to their mean values in redshift intervals of 0.5. For 
comparison, we also include the corresponding mean trends for 
the low metallicty sample (dashed line) and the complete galaxy 
population (solid line). 



Regarding the mean specific SFR (SSFR), defined as 
the ratio between the SFR and the stellar mass of a given 
galaxy (Fig. |4]), host galaxies are predicted to show SSFRs 
similar to those in the low metallicity sample, and higher 
than those in the complete galaxy population, for z < 2. 
Between z ^ 2 and « ~ 6, there are no significant differ- 
ences in mean SSFRs among the three samples, which have 
larger SSFRs that low redshift galaxies. This can be under- 
stood because, at higher redshift, galaxies have larger gas 
reservoirs which can feed stronger star formation activity 
and are, on average, less chemically enriched. At low red- 
shift, host galaxies seem to be particularly efficient at trans- 
forming gas into stars. Note, however, that the mean SSFR 
predicted by scenario II. 3 at low redshift is half an order of 
magnitude lower than the mean observed SSFR. 



4.3 Luminosity and colour 

Observations show that host galaxies tend to be bluer than 
the general population of galaxies observed at a given red- 
shift, and fainter than a typical L* galaxy. In fact, this trend 
can be nicely reproduced by our scenario II. 3, as shown in 
Fig. [5] The luminosities and colours of the predicted host 
galaxies are in excellent agreement with the observations 
compiled bySavaglio et al. (2009). From this figure, we can 
see that host galaxies are bluer than the complete galaxy 
population for z < 2, but have similar mean colours for 
higher redshfits. Our scenario predicts that host galaxies 
are more luminous systems in the _B-band, compared to the 
mean luminosity of the global galaxy population at all red- 
shifts. 



4.4 Metallicity 

In Fig. [6l we compare the mean cold gas metallicity of the 
host galaxies predicted by sce nario II.3 to di f ferent metal- 
licity estimations reported by ISavaglio et al.l ((20091). Since 



at z > 0.2 it is very difficult to distinguish HII regions with 
enough resolution, these authors generally measured the op- 
tical luminosity-weighted mean metallicity in a galaxy. For 
consistency, host galaxies at z < 0.2 were treated by them 
as the rest of the sample (i.e. integrating fiuxes over the 
whole galaxy). In this figure, we have incl uded metallici- 
ties o btained by using different indicators (see lSavaglio et al.l 
|2009| . for more details). In some cases, absorption lines in 
the optical afterglow can be used to obtain the metallicity 
of neutral cold gas along the line of sight of the GRB (the so- 
called GR B-DLAs). This is the case of 9 GRB-DLA systems 
studied by ISavagliq (|2006l ). all of them at z > 1.6. In this 
case, the metallicity could be associated more directly to the 
metallicity of the host galaxy, contrary to QSO-DLAs which 
are associated to HI clouds in the intergal actic medium. 
At z < 1, measurements of metallicity by ISavaglio et al.l 
(120091 ) are derived from hot gas and, in this case, the lower 
branch solution is preferred for the hosts. As shown in Fig.|Sl 
the predicted host galaxies metallicities are systematically 
higher than the GRB-DLA metallicities. This might indicate 
that LGRBs occur in regions of even lower metallicity than 
the mean me tallicity of t he co ld gas of the host galaxies. 
In this sense. iNuza et al.l (|2007l ). adopting similar hypothe- 
ses to generate synthetic LGRB populations, but using full 
cosmological simulations where the chemical enrichment of 
baryons was consistently followed with redshift, claimed a 
metallicity threshold of O.SZq for the progenitor stars in or- 
der to reproduce observations. More detailed information on 
the metallicity of individual stellar populations, which also 
takes into account the inhomogeneities of the interstellar 
medium, are needed to impro ve our understanding of this 
issue (see lPontzen et al.ll2009l and Art ale et al. in prepara - 
tion). Also along these lines, a work bv iModiaz et al.l (120081 ) 
compares the metallicity of GRB hosts associated to Type Ic 
supernovae (SN Ic) to broad-line SN Ic without GRB detec- 
tion at z < 0.25. The metallicity of SN-GRB host galaxies 
are estimated by computing the nebular oxygen abundance 
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Figure 5. Mean B-band luminosity (upper panel) and B — K 
colour (lower panel) as a function of redshift predicted by sce- 
nario II. 3 (dotted line) with its standard deviation (shaded grey 
band). Open d iamonds represent the observed host galaxies by 
ISavaglio et al.l (2009), while filled circles correspond to their mean 
values in redshift intervals of 0.5. For comparison, we also in- 
clude the corresponding mean trends for the low metallicty sam- 
ple (dashed line) and the complete galaxy population (solid line). 



by three different methods. In this case, they directly probe 
the environments of the progenitors because either they were 
measured at that location or they correspond to motallicities 
of homogeneous dwarf galaxies. Modjaz ct al. (200 ^) results 
are in agreement with the values obtained bv lSavaglio et al.l 
lJ2009f ) and with our own results at low redshifts. 



5 ENVIRONMENT OF LGRB HOSTS 

It is still under discussion if host galaxies are always lo- 
cated in regions of a given characteristic local density or 
if, dep ending on z, th ey inhabit different environments. Re- 
cently, [Cam£isi_et_alJ ([2009) estimated the cross-correlation 
function of their simulated host galaxies finding that these 
systems tend to map underdense regio ns. Observations 
provid e somehow contradictory results. iBornancini et al.l 
(^2004) f ound host galaxies to inhabit field regions. Lately, 
IWainwr ight et al.l (i2007) suggested that host galaxies might 
be biased towards interacting or merging galaxies. It is yet 
too soon for observations to provide a robust answer to these 
questions but theoretical predictions could positively con- 
tribute to this area. The advantage of our model is that 



Figure 6. Mean cold gas metallicity as a function of redshift 
predicted by scenario II. 3 (dotted line) with its standard deviation 
(shaded grey band) . Filled circles represent obse rved metallicities 
of host galaxies reported bv lSavaglio et al.l 1120091 ) , measured using 
the Te technique. We also include those measured by the R32 
upper branch (triangles) or lower branch (diamonds), and GRB- 
DLAs (asterisks). For comparison, we include the corresponding 
mean trends for the low metallicty sample (dashed line) and the 
complete galaxy population (solid line). 



taking into account the observability allows us to make a 
more robust prediction. 

In this section, we investigate the environment of the 
simulated host galaxies by means of three different estima- 
tors. The distance to the closest neighbour (di) of a host 
galaxy provides an estimation of the possibility of a host 
galaxy being interacting with another galaxy. The distance 
to the fifth neighbour (ds) of a host galaxy provides an es- 
timation of the density of the region the galaxy inhabits. 
Finally, the central halo virial mass (Mvir) is an estimator of 
the global potential well to which the host galaxy is bound. 

In the particular case of the closest neighbour, it is well 
known that p air interacti ons can enhance the star formation 
activity (e.g. iLambas et al. 2003 : Patton et al. 2005 ) and if 
LGRB progenitors are massive stars this would imply an in- 
crease in the probability of detecting a LGRB event. Since 
semi-analytical codes do not model tidally-induced star for- 
mation, we modified the SFR for galaxies with their clos- 
est neighbour located at distances smaller than O.lh^^ Mpc 
by a factor of t wo, motivated by observational and numeri- 



cal resuhs (e.g. iLambas et al.ll2003l : IPerez et al.ll2006l . l2009l : 



iDi Matteo et al.l 120081 ). However, this correction does not 
change any of the trends reported in this work. But the role 
of interactions should be further tested with models that 
properly include their effects. 

We estimate the cumulative distribution of these three 
environmental estimators (di,d5,Mvir) predicted by sce- 
nario II.3. To quantify behaviours as a function of redshift, 
we determine the values di,5o, ds.so and Mvir, so at which the 
corresponding cumulative fractions reach 50 per cent in the 
case of di and ds, and 80 per cent in the case of Mvir. In 
Fig. [71 (upper panel), we show di,5o as a function of redshift 
according to scenario II.3, together with the trend for the 
low-metallicity sample and the complete galaxy population. 
At high redshift, the three samples tend to have the near- 
est neighbour at approximately similar distances. We only 
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note a slight trend for host galaxies to systematically have a 
closer first neighbour than galaxies in the other two samples 
for z > 1 . From z ~ 1, low metallicity galaxies tend to have 
their closest neighbour further away than either a galaxy in 
the complete population or a host galaxy. The latter tend 
to inhabit lower density regions than galaxies in the com- 
plete population. Our models indicate that observable host 
galaxies would only have a slightly higher probability to be 
in a galaxy pair compared with the general galaxy popula- 
tion at high redshift (i.e. 50 per cent of galaxies have the 
first neighbour closer than ~ 50kpc). 

The global environment is quantified by ds.so. As shown 
in Fig. [7| (middle panel) at high redshift (z > 2) , galaxies in 
the three samples reside in similar environments. For lower 
redshift there is a clear trend for galaxies in the low metallic- 
ity sample and host galaxies to inhabit lower density regions 
than galaxies in the complete galax y population. This re - 
sults agree with the numerical work of lCarnpisi et al.l (120091). 
and al so with the observational findings of lBornancini et al.l 
(I2OOJ), taking into account that the mean redshift of LGRBs 
observed at that time was (2) < 2. 

In the case of the virial mass, as shown in Fig. [7] (lower 
panel), we find that 80 per cent of observable host galaxies 
have haloes less massive than iq^^o-^^-^ ^^ g^ g^j^y j.g(j_ 
shift. As the structure forms and groups and clusters ag- 
gregate hierarchically, galaxies tend to inhabit larger dark 
matter haloes. However, host galaxies stay within a nar- 
rower range of mass haloes. The inversion in the relation 
observed in this figure can be understood considering that 
host galaxies observability depends strongly on the star for- 
mation a ctivity which, in turn , depends strongly on envi- 
ronment (jPoggianti et al.l 12009 ). Hence, from z ~ 2 active 
star forming galaxies (i.e. those which are important gen- 
erators of LGRBs) with mean stellar masses of ~ 10^° Mq 
reside again in slightly smaller dark matter haloes (Fig. [2] 
and Fig. [3}. Note that this is not the case for the general 
low metallicity sample which although residing, on average, 
in smaller haloes than the global galaxy population, tend 
to systematically inhabit larger ones with decreasing red- 
shift, as expected in a hierarchical scenario. The observabil- 
ity condition which is closely linked to the star formation 
activity produces this kind of halo downsizing scenario for 
host galaxies at low redshift, while the mean stellar mass 
remains approximately constant. 
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6 CONCLUSIONS 

Our results indicate that the observed LGRB host galaxies 
properties can be reproduced assuming that LGRB progen- 
itors are massive stars and occur in galaxies with moder- 
ately low mean cold gas metallicities. The joint requirement 
for the synthetic LGRBs to reproduce the observations of 
BATSE and for the simulated systems that host them to 
reproduce the stellar mass distribution of the observed host 
galajcies determines the minimum mass for the progenitor 
star (m > TSMqJj and the maximum metallicity cut-off 



■^ The m ass thres ho ld fo r LGRB progenitors was computed 
for the [C habrier (2003) IMF used in the cata logue of 
iDe Lucia fc Blaizot (2003), and would be 54 Mq if a ISalpeten 
(ligSSl ) IMF had been used, instead. 



Figure 7. Environmental study of host galaxies predicted by 
scenario II. 3 (diamonds), the complete galaxy population (stars), 
and the low metallicity galaxy sample (crosses) . Upper panel: Me- 
dian distance to the closest neighbour (di_5o). Middle panel: Me- 
dian distance to the fifth neighbour (^5,50)- Left panel: Virial 
mass of the central halo at the 80 percentile (Mvir,8o)- 



(0.6 Z0) for the cold gas metallicity of the simulated host 
galaxies. Our scenario II. 3, which satisfies these conditions, 
succeeds at reproducing the dependence on redshift of stel- 
lar mass, luminosity, colour, SFR, SS FR and metallicity o f 
the observed host galaxies compiled bv lSavaglio et al.l (|2009l ) 
over the redshift range < 2 < 3. 
Our main findings are: 

(i) The average stellar mass of host galaxies is higher than 
the average stellar mass of the complete galaxy population. 
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and remains within a narrow mass range around 10^*^ M© , 
with a very weak trend to lower mass for higher redshifts, 
due to the double requirement of low mean gas metallicity 
and high star formation rates. We note that the dispersion 
is approximately 0.5 dex around the log mean stellar mass. 

(ii) Compared to the characteristic st ellar mass M* ^ 
-^qIO.63 y^^ estimated bv lBell et aP l|2003h . host galaxies are 
low stellar mass systems (88 per cent of our host galax- 
ies have stellar masses between ~ ]^qS-5-io.3 ]y[g) However, 
compared to the mean stellar mass of the complete galaxy 
population, host galaxies tend to be massive galaxies. 

(iii) The SFR of host galaxies is larger than that of galax- 
ies in the complete population at all redshifts. At z > 2, 
the SSFR of observable host galaxies is comparable to the 
complete galaxy population, but at lower redshift it is sys- 
tematically higher. As a consequence, host galaxies are, on 
average, bluer than the global galax;y population. The dis- 
persion found for the properties of the observable host galax- 
ies reflect the different histories of formation of galaxies at 
a given stellar mass. 

(iv) Our results support the claims for a metallicity 
threshold to reproduce the properties of the current observed 
host galaxies. The global metallicity threshold of 0.6 Z© we 
derived from our models is an upper limit for the metallicity 
of the LGRB progenitor stars in the case of no dust effects. 
The comparison of the mean metallicity of host galaxies 
predicted by our model with GRB-DLAs observations sug- 
gests that LGRBs might be produced in stars o f even lower 
meta l licity as pointed out by other authors (e.g. lNuza et al.l 
120071 : ICampisi et al.l 120091 ). At low redshift, our results are 
in agreement with S N-GRB local metallicity estimations by 
iModiaz et al.l (|2008r ). However, observations show a large 
spread in metallicity and there are also observational un- 
certainties that can affect its determination. Hence, results 
should be taken with caution. 

(v) For 1 < 2 < 6, host galaxies seem to be slightly more 
likely to be in pairs than galaxies belonging to the other 
samples. There is, however, a change of behaviour for z < 1, 
where the at least 50 per cent of the host galaxies seems 
to have their closest neighbour further away than galaxies 
in the complete population, but closer than galaxies in the 
low metallicity sample. We highlight that we are considering 
neighbouring galaxies of a mass above 10^ Mq . This limit is 
set by the numerical resolution of the Millennium Simula- 
tion. Lower mass companions could also imprint morpholog- 
ical perturbations and trigger star formation activity (e.g. 
iLambas et al.ll2003l ) but we cannot pursue this analysis fur- 
ther on with this galaxy catalogue. 

(vi) Regarding global environment, our model suggests 
that, at z > 1, observable host galaxies would preferentially 
inhabit environments of similar density to those populated 
by the general population and of slightly higher density than 
those inhabited by the low metallicity sample. Towards z — 
0, observable host galaxies tend to be progressively located 
in less dense environments, which becomes as subdense as 
the regions where l ow metallicity galaxies reside (see also 
ICampisi et"aLll2009l ). 

(vii) Our results suggest that observable host galaxies 
tend to have dark matter haloes in the range iq^^o-^^-5 Mq, 
regardless of redshift, and show a slight signal for halo down- 
sizing from z ~ 1, distinguishing them from the other two 
samples which follow the expected halo mass growth in a hi- 



erarchical scenario. This result is consistent with observable 
host galaxies being systems with mean stellar masses ap- 
proximately constant, regardless of the age of the Universe. 

Our results are mainly a consequence of the joint re- 
quirements to have high star formation activity to ensure 
observability and to reproduce the current distribution of 
stellar masses of observed host galaxies. However, within the 
current constrains, galaxies with high masses would be too 
metal-rich to produce LGRBs and low mass systems would 
have low probability of being observed. If the observed sam- 
ple were modified by the incorporation of other galaxies, 
for example dusty hosts, then our model would need to be 
readjusted accordingly. 
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